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Abstract 
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b se^e~, b s'^ViVi in supersymmetric model with bilinear R-parity 
violation. We compare our calculations with the evaluations in SM and 
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I. INTRODUCTION 


It is being increasingly realized by those engaged in the search for the supersymmetry 
(SUSY) that the principle of R-parity conservation, assumed to be sacrosanct in the 
prevalent search strategies, is not inviolable in practice. The R-parity of a particle is dehned 
as R = (—i)^+3S+2S^ pgg Pg violated if either baryon (B) or lepton (L) number is 
not conserved in nature, a fact perfectly compatible with the non observation of proton 
decay. Under R-parity violation the phenomenology changes considerably the most 

important consequence is that the lightest supersymmetric particle (LSP) can decay now. 
However, the way in which R-parity can be violated is not unique. Different types of R-parity 
violating interaction terms can be written down, leading to different observable predictions. 
In addition, R-parity can be violated spontaneously, instead of explicitly, whence another 
class of interesting effects are expected IQ. If the phenomenology of R-parity breaking has 
to be understood, and the consequent modihcations in the current search strategies have to 
be effectively implemented, then it is quite important to explore the full implication of each 
possible R-breaking scheme. 

The R-conserving part of the minimal supersymmetric standard model (MSSM) is spec- 
ihed by the superpotential 

y^MSSM = + EijdijHlQ^jD'^ 

+e,,ujjHfQp-^, ( 1 ) 

where /, J = 1, 2, 3 are generation indices; i, j = 1, 2 are SU(2) indices; and e is a 
completely antisymmetric 2x2 matrix, with £12 = 1. The ”hat” symbol over each letter 
indicates a superheld, , H^, and being SU(2) doublets with hyper-charges — 1, 

— 1, and 1 respectively; U, D and R being SU(2) singlets with hyper-charges — | and 2 

respectively. The couplings u/j, du and Rj are 3x3 Yukawa matrices, and fi is parameter 
with unit of mass. If now the bilinear R-parity violating interactions are incorporated, the 
superpotential takes the form P 
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W = Wmssm + VVl 


(2) 


with Wl = and e'j is the parameter with unit of mass. 

The soft SUSY-breaking terms: 

Csoft = - mirL]*Ll 

- mliQl*Ql - 

-mljU^*U^ + (miAfiAs + msA^A^ 

RttisXqXq + h.c.) + (^B fiEijR^ + BjejSijHfLj 

RSijlsi^HlVjR^ + EijdsifiHlQjD^ 

+EijUsifiHiQjU^ + h.c.) (3) 

where rrij^i, mjj 2 , rri^i, ttIqi, mj^i and rrifji stand for the mass squared of the scalar 
fields while mi, m 2 , m 3 denote the masses of the SU(3) x SU(2) x U(l) gauginos A^, A^ and 
Xb, B and Bj (/ = 1, 2, 3) are free parameters with unit of mass. 

Wl and lepton-number breaking terms in Cgoft give a viable agent for R-parity breaking. 
It is particularly interesting for the fact that it can trigger a mixing between neutralinos and 
neutrinos as well as between charginos and charged leptons, resulting in observable effects. 
For 63 can reach a large value (such as jcgl ~ SOOGeV) even though we have considered 
the experimental fact that rrir = 1.77GeV and m^,^ < 24MeV |^, the mixing between 
charginos and charged leptons may give an important contribution to the rare B decays 
such as 6 —>• s -1- 7 , 6 —*• se’''e“, b ^ 

On the other hand, flavor-changing neutral current (FGNG) processes have constantly 
played a major role in the development of electroweak theories. FGNG studies have found 
kaon physics as their best research ground so far. However, the achievements in B-physics 
make B-mesons the new challenging frontier in the study of FGNG phenomena. At least 
some of the rare B-processes which are induced by the FGNG transitions b —>■ s and b —>■ d 
are within the reach of present machines. Needless to say, the test of such rare B-processes 
presents a probe of the validity of crucial ingredients of the Standard Model (SM) and. 
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possibly of the existence of new physics beyond the SM at the low energy scale. 

b ^ s + j and b —^ se~^e~ in the SM are analyzed in Ref The SM prediction for 

the branching ratio (Br) of the inclusive decay 6 —*• s + 7 , once the large QCD corrections 
are included 0, is a few times 10“^. As for the semileptonic b se~^e~ decay, the QCD- 
corrected SM prediction for its Br is about 10“® [^]. 

FCNC rare B-processes have been widely analyzed as a potential probe for extensions 
of the SM implying new physics at the TeV scale. In the model with Two Higgs Doublets 
(THDM) and no tree level FCNC, the rare B decays and Bd — Bd mixing are computed in 
Ref P, ^ . As for left-right symmetric models, the simplest choice of taking the CKM mixing 
angles in the right-handed sector equal to the corresponding left-handed ones (manifest or 
pseudo-manifest left-right symmetry) does not give any appreciable effects with respect to 
the SM estimates of rare B-processes Analogous results are obtained in the general 

case if no hne tunings of the parameters are allowed . 


FCNC in the supersymmetry with R-parity has been discussed widely in Ref |]^, Ref 
, FCNC in supersymmetry without R-parity has been discussed in Ref []^, the review 
of FCNC in the supersymmetric models can be found in Ref [^. 

This paper is organized as follows. In sect. II, we give a description of the general 
structure of the supersymmetry with bilinear R-parity violation (BRPV). In sect. Ill, we 
provide the complete analyses to 6 —> s -f 7 , 6 —*• se~^e~, b —> sY^ViVi in supersymmetric 
model with bilinear R-parity violation. In sect. IV, we perform a numerical analysis of those 
processes and compare them with the SM predictions and the experimental results. For 
completeness, we give the analytic expression for Feynman integrals which one encounters 
in the evaluation of the amplitudes listed in sect. Ill in appendix. A. A systematic notation 
for the relevant Feynman rules involving R-parity breaking terms are introduced in appendix. 
B. In appendix. C, we give the mass matrix of r neutrino-neutralinos. For simplicity, we 
neglect the generation mixing of slepton and squark in the following analyses. 
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II. MINIMAL SUSY MODEL WITH BILINEAR R-PARITY VIOLATION 


From the superpotential Eq. (|^), we can perform an operation that is the same as in the 
standard model by redefinition of the fields 


vb^bb 

Ll Vi;Li, 

Fb (4) 

One can diagonalize the matrices Ijj, ujj and djj, the superpotential has the form 

W = + hSijHmR^ - 

-HlQ{)U^ - di{HlQl - HlC^JQi)b^ 

(5) 


where C is the Kobayashi-Maskawa matrix and e/ (J = 1, 2, 3) have the definition as: 

c = r^VQ.. 

£; = «'W"- (6) 


The soft breaking sector has the form 

- mliV*U - m%R^*lb 

—rriQiQbQl — m\,ib^*b^ — rnbiU^*U^ + (miAsAs 
+7712 A^A^ + tti^XqXq + h.c.) + + BjejEijH^Lj 

+eibsitiHmR^ + dsifx{-HlQi + C^^HlQb)b^ 

+ HlQ{)U^ + h.c.}. (7) 


For simplicity we take from now on ei = 62 = 0, in this way, that only r-lepton number is 
violated. The electroweak symmetry is broken when the two Higgs doublets and the 

r-slepton acquire vacuum expectation values (VEVs) 
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(8) 


i_ + vi + 

[ m j 




(9) 


( 10 ) 

V y 

Note that the gauge bosons W and Zq acquire masses given by rri^ = \g‘^v‘^ and m| = 
1(5'^ + where + n| and g, g' are coupling constants of SU(2) and U(l). 

We introduce the following notation in spherical coordinates 


vi = V cos6*„ cos /3, 

V 2 = V sin l3, 

V 3 = V sm9^ cos p. (11) 


When the angle 9^ equals to zero, this sector will change back to the MSSM limit exactly. 
The massless neutral Goldstone boson can be written as: 


= cos 9^ cos (3^Pi — sin /3(p2 + sin 9^ cos (12) 

In the model with bilinear R-parity violation, the charged Higgs bosons mix with the left 
and right handed r-slepton. In the original basis, where <hc = {H 2 *, Hf,fl,fR), the scalar 
potential contains the following mass term 

eg = (13) 

with the symmetric matrix is given by (here the matrix is too big to be written in full 
so we write it by each element individually) 
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cl,2 
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cl,3 
2 

cl,4 
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c2,2 
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c2,3 
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c2,4 

2 

c3,3 

2 

c3,4 

2 

c4,4 


y / 2 


<^3 I rj "^^2 

€ 3/1 -h BfX —, 

Vi Vi 


2 \ I 2/2 I ^3 

■^ 3 ) + 2'^3^3 + *'3/^ — 

+ 5/i, 

1 / 

-tiitis - e^iJ, - -I3V1V3, 

V2 . , /it^3 

jK + ^^ 1 ) - ^ 363 — + 5/i—, 

4 112 ^2 


-112^^3 - ^ 363 , 

|/^^3 + ^eatii, 

/ 2 2 n I 

-K -'^^l) +^3- 

V 3 


R e3t^2 /i 2 

■ ^3-K 17 ^^ 1 , 

113 2 


= /3/it;2 - 


72 




.^,2 


This matrix has an eigenstate: 

G+ = ^ Zlfil 

i=l 

= -(^^ 1 ^ 2 * - V 2 HI + V3fl) 

V 

= COS 6^ cos PHl* — sin l3Hl + sin cos jSf^ 


(14) 


(15) 


with zero eigenvalue, it is the massless charged Goldstone boson. In the physical (unitary) 
gauge, are absorbed by bosons and disappear from the Lagrangian. The other three 
eigenstates , fi, f 2 can be expressed as: 

H* = Y.Z]tK, 

i=\ 

n = i: 4 T. 

i=l 

= ( 16 ) 

i=l 

Similarly to the Higgs bosons, charginos mix with r lepton forming a set of three charged 
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fermions r, ^2 > ^3 [ 0 ) 0 - Iii fhe original basis where = (—iA"*", H^, t^) and pj ^ 
Hf, rh), the charged fermion mass terms in the Lagrangian are 


= —i’ ^M.fljj^ + h.c. 
with the mass matrix is given by |T3|JT^: 


(17) 



2m2 

ev2 

\/2s\y 

0 ^ 

Mf = 

ev-\ 

\/2sw 

/i 

hv?, 

U2 


eu3 

\ \/2sw 

£3 

I3V1 

V2 / 


(18) 


where Sw = sin^w and ^ . Also the two mixing matrices and Z~ appear in 

the Lagrangian, they are dehned by the condition that the product {Z^)"^M.fZ~ should be 
diagonal 


{Z^fMfZ- = 


^ m-- 0 0 ^ 

0 m-- 0 

^2 

0 0 m-;- 


(19) 


Here, we assume m-- 


rur and m~- > m~- > m-- 

^3 ^2 


III. COMPLETE ANALYSIS OF THE RARE B-PROCESSES IN SUSY MODEL 
WITH BILINEAR R-PARITY VIOLATION 

Because we neglect the generation mixing of sleptons and squarks, there are no the 
couplings such as 7^ J)^ ^hf^PL,RqdDl 2 (/ ^ J), R^'yf.P^Re^{I 7^ 

J), XcIt^ELMiP 1,2 {I 7^ J) and \Q'y^PL,RqdDi 2 (-^ 7^ J) where /, J = 1, 2, 3 are the 
generational indices and k^, k^, Xq represent the neutralinos and gluinos. Under 

this assumption, the contribution of gluino and neutralino is zero when we compute the rare 
processes such a.sb—>s + 'j,b^s + e+e“ and b ^ s + J2 Vi^i at one loop level. 

In this section, we will use the effective Hamiltonian theory to discuss those processes. 


The method of the effective Hamiltonian theory was hrst used by Ref [0 and has been 
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developed over the last years ||2^ [Q. It is a two step program, starting with an operator 
product expansion (OPE) and performing a renormalization group equation (RGE) analysis 
afterwards. The derivation starts as follows: if the kinematics of the decay are of the kind 
that the masses of the internal particles rrii are much larger than the external momentum 
p: then the heavy particles can be integrated out. This concept takes a concrete 

form with the functional integral formalism. It means that the heavy particles are removed 
as dynamical degrees of freedom from the theory, hence their helds do not appear in the 
(effective) Lagrangian anymore. Their residual effect lies in the generated effective vertices. 
In this way an effective low energy theory can be constructed from the full theory. In the 
framework of the standard model, strong interactions are known to give sizable contributions 
to FCNC processes, the inclusion of the QCD corrections increases the electroweak rate by 
about a factor two for 6 —sy process, and enhances the rate about 30% for b se'^e~ 
transition [^]. Now, let us derive the effective Lagrangian for b —> sy, b —> se~e'^ and b —> 
sS GG first (at the scale), then we analyze the renormalization group equation (RGE) 
in the leading logarithmic approximation (LLA), i.e. to sum up the terms ln('^^)]’^ to 
all orders n (n=0, 1, ■ ■ cx)) in perturbation theory. 


A. The bsZ, hs^ couplings and box diagrams in SUSY model with bilinear R-parity 

violation 


We follow the method that was used in Ref []^. In the't Hooft Feynman gauge, the 
one-loop diagrams for the induced bsZ coupling are shown in Fig. ^ The diagrams Fig. 
|^(a). Fig. Q(b) and Fig. |^(c) represent the SM contribution of the 6 —>• s transition. In each 
group, the first two diagrams (except for Fig. |I|(c)) are self-energy part. Since the weak 
current is not conserved, we only need to expect a non-vanishing zeroth-order contribution 
in the momentum q (q denotes the momentum transfer through the Z-boson). The induced 
bsZ coupling takes the form: 




(20) 


9 








where i = a, b, c, d, e and Pl,r = The F* can be written as (at the scale): 


r'“’ = „ . L - ^c;,c,i{d -ism"«„,)| 2 /f (i,„)- 

[4:71)^ sin^ Ow cos Ow 2 3 

4^f3b\xtw)] - -^sm‘^ew)f3bixtw) - 

2 3 

- sin^ OwXtwfSb\pw) + - cos^ OwffJixtw)}, 

(47r)^sin c/y/cosc/y/ 2 3 2 


(21) 


^tw 


^(cos^ 9w - sin^ 9w) {xt^)}, 


ph) = 


pW = 


(47r)2 sin^ 9w cos 9w 


C,’,CttV2 siii'i 0wXtu,f!4’ ixtv,), 


( 22 ) 

(23) 


(47r)2 sin'* 9w cos 9w 


C;.C,i.{(i - hin"9».) 


^tw 


Zz];zf[fP(x,.,)- 


‘^fslHxtSi)] + 77-7 


Xtw 


-2 3 

4 


2 sin 9^ sin (5 ^ 


2 sin 9^ sin f3 ^ 

isiri^«„,4°hi..)l - 




4 4 


Y.T.ZhZT\{^ox‘0„,- 


P(e) = 


8 sin 9^ sin /3 ^ ^ 

i=ij=i V2Mwsm9ySmp 

2 3 


(24) 


(47r)^ sin^ 9w cos 9w 


'72i* '7+ 
^Q3 ^ 2 i 


)r(;- 


1 inl :4 2 i rntZ^^^ZP 

-s..^.„,)l/f (.,-,) - 2 /i;p,-,)l + ZU-z$zi + 

V+* I ,, ^ ryll* r;xli ^ ■ 2n f r-- - ) I 




3^1 , mtZgZ^* 

EEk-^ozZ,r+ 

i,j=i 1=1 ^ 




\/2Mw sin 9y sin (3 
1 
2 




+ 25iiCos20H2)^/iP(x--,-,x--,-) +x--,-x--,-(Z+*Z+ 


\/2Mw sin 9y sin /3' 

+* Y + 


+2Sij cos 29 w)fS} {x-^-i, x--r)]}, 


i * 1 ’ 


(25) 


2 

where Si=i, 2 , 3,4 = {G~, H~, fi, f 2 ) and Xap = Zq 3 is the top-squark mixing matrix, its 
dehnition can be found in appendix. B. 

Notice that the ultraviolet divergences cancel separately for each of those equation. In 
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obtaining the form shown, the nnitary property of Zh, Zq3, and Z has been nsed, 
together with the nnitarity of the Kobayashi-Maskawa matrix. 

The compntation of the photon exchanged contribntion is somewhat more involved and 
reqnires the calcnlation of the indnced bs'j conpling np to second order in the external 
momentnm. The diagrams need to be compnted are those of Fig. with Z being replaced 
by 7 . The indnced bsj conpling takes the form (at the Mw scale): 

r« = - qj)p^ + F^^-f^{msPL + (26) 


where i = a, b, c, d, e and 




(dvr) 2 sin^ 6w 


^tb^ts 






= 


^ ™ 1 .2 70 ) 


(dvr) 2 sin^ 6^ 


CtbCt 


2 a 




w 


^f5d\^tw) + ^f3a\xt^) - - 


In 


Ml 




mt 


Pf ^ = 
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Pf^ = 


(dvr) 2 sin^ 9w 

(d7r)"'sm Mui 2 


3a; 


tw n(l) 


''Ml 

d 2 ). 


12 


~^ffdi.Xtw) + 


Ub i^tw) - -^fL ixtw) - 


12 


fS\xtw)}, 


= 


= 




(d7r)2 sin^ 6w 




'Ml 


(dvr) 2 sin^ 9^ 

e3 ^ 


-CtbCl 


dM^ 


fil\xtw), 


(dvr) 2 sin^ 9w 




^721 '72i^ 




7 6M^. siii^ Ow siii^ /3 ^ 3 


(27) 


(28) 


(29) 


(30) 

(31) 

(32) 
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(33) 




pid) _ 
^2 — 


1 


(47r)2 sin^ 6w 

'72i '72i* 

/( 2 ) 


Sin 




a:f.. 


^ 6M^ sin^ 6 *vk '■' sin^ {3 

Xts, 

w 


' 72 * ' 72 ** ' 71 ** ' 72 * 


■e 


sin (3 cos (3 


)/iP(a;tsJ 


' 72 * ' 72 ** 

Al)i 




Z li* ^2i 
H 

2 sin [3 cos f3 


' 72 * ' 72 ** 


3sin"/3' 


' 72 * ' 72 ** 


6 sin (3' 


(34) 


= 


mtZ^3^Z+; 


/O \ ^ \ ^ ■*■ I r^li /T’H-* 


(47r)2 sin 


*=1,2 i=l 


y/2Mw sin 6y sin {3 9 


nh'i'A-i.) + 


(35) 




Fi"^ = 


(47r)2 sin 6w 


0*c;, i: 


/Ol/f • *3 • 9'^4 c (^k. 1*) 

v2Mvi/sin6't,sinp 2 j 


*=1,2 i=l 








:) 




\/2Mw sin 9v sin (3 ^/2Mw sin 6^ cos (3 




c(l)/ 


(36) 


The function f(s) is defined as 


_ _2 _ 2 ^ I S if s < z, 

3 F|(i + i)yr3T|i„i±^_i^], if 

where and /(O) = 1. The other functions are given in appendix. A. 


(37) 


The box-diagrams that contribute to the b —>■ se’''e are shown in Fig. |^. The effective 
Lagrangian takes the form (at the Mvk scale): 


A^^^isYPMej.PLe) 


(38) 


with 


^box _ 


(47r)^ 4sin^ 9w 


/O* /O 

^ts^tb 


Mh 


/«(*"«. 0). 


(39) 


J^OX 


(47r)2 4sim 9w 


3 2 1 

n* n _' 7 +* ' 7 +(' 7^^* 7 + 




•17 1 

2J = 1 / = 1 ti 


\/2Mw sin 9^ sin [3 
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(40) 


m,zl,z+c 


[-Z^- z^* + 

\/2Mw sin 6^ sin j3 ' 




As for b —>■ sUeiZe, it is analogous with the case of 6 —*• se~^e~. The box diagrams that 
contribute to the b —>■ sugize are given in Fig. The effective Lagrangian takes the form (at 
the yiw scale): 


Bt°^{srPLb){9gj^PLUg 


(41) 


with 


Qbox _ 


^box _ 


1 ,(i). 


(47r)2 4sin^ 

1 p4 1 3 

^ /O* /O \ ^ \ ^ \ ^ rylh fylh* ry— ry — ^ 


(42) 


m,z^A;z+ 


_ yV-* _j_ 

y/2Mw sin 6y sin (3 


|2r(i) 


/ie (%-iV%-tva:g-4-) 


(43) 


where Z^i is the mixing matrix of slepton (1=1, 2 is the index of generation), its dehnition 
can be found in appendix. B. 

B. The width of 6 ^ s + 7 in SUSY model with bilinear R-parity violation 

The total amplitude of the decay 6 —»• s + 7 can therefore be written as: 


AUb ^ s + ^) = F.OIj, 


(44) 


with F 2 is the sum of Eq. (p^), Eq. (|30|) , Eq. (|32[) , Eq. (|34[) , Eq. (^). Where O'ln = 
and the contribution of (9^^ = mge^s^'y^PLb is neglected since it is order of 

O(^). 

By denoting the total amplitude at a scale p as E2(p), the QCD- corrected amplitude at 
the scale of the process (~ rrib) is then given by []^ 

n(m,) = r^‘{F,(Mw) + - 1) + - 1)]}. (45) 

135 189 


where 
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1 


(46) 


7] = as{mb)/as{Mw) ~ 1 - 8 , 


— 

^2 — 


-CtbCl 




{4:71^ 2sm^ 

The property C*gCcb ~ —C^gCtb for the 3x3 CKM matrix has been used in the previous 
equation. The inclusive width for 6 —> s + 7 decay is hnally given by 


ml 


T{b^s + -f) = ^\F,{mb)\ 


(47) 


Where we have neglected the phase-space factor of order O(^). We calculate the corre¬ 
sponding branching ratio as in Ref by making use of the semileptonic decay b —> ceu, 
one gets: 


Br{b —>• s 4- 7) = ^ Br{b —*• ceu) 


(48) 


r(6 —>• ceu) 

where for Br{b —^ ceP) we use the averaged experimental value 0.11 [^]. The QCD 
corrected width for the semileptonic decay b —> ceP is H 


Tib^ceP) = ^p(^,0,0)|a7{l - r:77£i/(7£,0,0)} 


2 as{mb) 

dvr ^mb 


(49) 


where the phase-space factor p is 0.447 and /(^, 0, 0) = 2.41, Gp is the Fermi constant. 


C. The width of 6 —> se'’'e in SUSY model with bilinear R-parity violation 

This decay has been often considered the benchmark of charmless b-decays with strange 
particles in the final state. We will provide below the amplitude for b —> se+e“ including 
QCD effects. In what follows, the same conventions as in the previous decay have been 
adopted. 

We begin by considering the diagrams which induce the effective flavor-changing coupling 
of the photon to quarks (photon penguins). They are given by the diagrams shown in Fig. 
with a lepton line attached to the photon propagator. We consider separately the monopole 
(LL) and dipole (LR) form factor, which are related to different effective operators. 
(i)Photon penguins (LL component, at the scale) 
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•AJot^ih se+e ) = cFiOllv 
= F^Ollv 


(50) 


where Ollv = {s'y^PLb){e'yf^e) and Fi is the sum of Eq. (pTl) , Eq. (^, Eq. (|3lD , Eq. 
Eq. (HD. 

(ii)Photon penguins (LR component, at the scale) 




se+e ) = eF 20 LRv 
= F2 Olrv 


(51) 


with F 2 is the sum of Eq. (1^) , Eq. (H)) Eq. (|3^ , Eq. (H), Eq. (HD and the operator Olrv 
is dehned as: 

Olrv = mb^{s4^^PRb){e'y^e) 

(iii)Z°- penguins. The process b se~^e~ is also induced by the effective EC coupling of 
the to quarks. The total amplitude coming from the -penguins can be expressed as 
(at the Mvy scale) 


Afotib se+e ) = 


Mz sin 9w cos 9w 


^z{—-Olll + sin^ 9wOllv) 


= A (--OlllF sin 9 wOllv) 


(52) 


where the new operator Olll is given as 

Olll = {srPLb){e'y^PLe) 

and Tz = Y1 is the sum of Eq. (|2T|), Eq. (|2^),Eq. (HD) Eq. 

(iv)Box diagrams. The relevant contributions are given in Eq. 


and Eq 
, Eq 


In order to implement the QCD corrections, let us rewrite the total amplitude at the 
Mw scale as: 


nb^se+e ) 
•A~\- 


Hot 


{Mw) — Allv{Mw)Ollv + Alrv{Mw)0 lrv + Alll{Mw)Olll (53) 
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with 


Allv{M\y) — Fi + siv? ^ 

Alrv{Mw) = F^, 

Alll{Mw) = -A ^/2 + ( 4 - + 


(54) 


Renormalization at the mb-scale leads to [22 


ALLvimb) = A1 ^y{Mw) + sin^ 9wA^{Mw) + A!q ^ ~ 

-|[1 - T®)} - + fMP’f* - 

#1. 

. / N -!§. r A / 7. ^ \ wrlie ,10 ^ . 58 

ALRvi^b) = V ^^{AiRviMw) + “ 1) + ~ 

ALLii.'mb) = Alll{Mw), 


29, 

23 


(55) 

(56) 

(57) 


where Aq = F^^/Fiia and r/, F 2 have been dehned in Eq. 
The differential decay rate is then given by ||15||: 


+(1 + -)\Alrv\'^ — ^Re[{2ALLv + ALLL)A*i^j^y]} 
s 


where s = -^ 

mf 


( 68 ) 


D. The width of 6 —> in SUSY model with bilinear R-parity violation 

The transition b ^ sJ2 is induced by Z°-penguins and box diagrams, which, at the 
leading order, lead to the same effective operator. The peculiarity of this process is that it 
is not affected by the QCD renormalization. This is simply understood by noticing that the 
current {sj^Pib) is conserved in the limit of vanishing quark masses. Conserved currents 
have canonical dimensions and no divergent counterterms arise. This ultraviolet behavior is 


not spoiled by consideration of hnite quark masses . 
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The original electroweak sensitivity to the top mass is therefore preserved. On the other 
hand, the experimental search for this rare B- transition is understandably much harder 
than for the previous semileptonic decay. 

For the e-neutrino and //-neutrino, we obtain the following result 


Atot{b suiUi{i = e, //)) = (59) 

where A'^ = + Bl°^ with A^ = . „ —5—F^ and Bl°? is given in Eq. (0^) and 

' a ' 0 smOyycosOyy ^ a,b o 

Eq. (^31). The decay rate for the e, //-neutrino is then given by 

c; 

777 ^ 

T{b suiPi{i = e,//)) = ^535^3 (60) 

As for the r-neutrino, the case is complicated because the r-neutrino mixed with the 
neutralinos under our assumption. The interaction between the Zq boson and r-neutrino 
can be written as: 


^ Zq — h'r — h'r 


-Z^Ur^^PLl'r{\Zji 


3,1|2 


74,1|2 

\ 


+ I4T}. 


(61) 


2 sin 9w cos 9w 

where Ztv is the mixing matrix of r-neutrino and neutralinos, its definition can be found in 
the appendix. C. For the r-neutrino, we have 


Atotip >■ SUrl^r) — FmixA 


(62) 


with 


J- rr 


173,112 174,112 I 175,112 

\^N I \^N I + \^N I • 


(63) 


The decay rate for the r-neutrino is given as: 


F(6 —i> SUr^r) = 


ml 




(64) 


153671^ 

here we neglect the mass of r-neutrino. We may now sum over all three types neutrino 
species: 


F(6 ^ s^ViVi) = 


ml 


15367/3 


{2 + 


(65) 


As we shall see in the next section this leads to a rate for b ^ sJ2 which is about four 
times larger than the rate for b —»• se'''e“. 
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IV. NUMERICAL RESULTS 


In this section, we will compnte the branch ratio of the processes that have been analyzed 
in sect. III. In the numerical evaluation below, we take a = ~ T 2 IT’ ^s{Mw) ~ 0.118, 

rrie ^ O.SllMeV, 171 ^ ^ 1.777GeV, Mz ~ 91.187GeV,Mvi/ 80.330GeV, rrit ~ 174.GeV, 

rric ~ 1.3GeV and mb ~ 4.3GeV. 

In order to hnd out the allowed region in the parameter space, one has to take a number 
of constraint into account. First, we note that mi, m 2 , 63 , p, tan/?, tan 6 '^ and I 3 are 
the parameters that enter into the chargino and neutralino mass matrices. The strongest 
constraint on them follows from the fact that the r mass has been experimentally measured 
therefore, for any combination of those parameters, the lowest eigenvalue of Eq. 


should agree with m,-. Also, has a laboratory upper limit of 24MeV on its mass. The 
two restrictions, together with the positive-dehnite condition of the mass squared matrices, 
constrain the parameter space in a severe manner. 

We open our discussion by considering the experimental results of rare B processes impact 
on the mass spectrum in the model. At present, the experimental bound on the b s'j, 
b —>• are 2 x 10“"^ < Br{b —>5 + 7 ) < 4.5 x 10“"^ and Br{b —> < 

5.7 X 10“® [Q respectively. Not having the generality lost, we assume = I 3 in the 
numerical calculation and the value of can be determined from Det\m^ — AdJ-Al/l = 0 
when the relevant parameters are given. Furthermore, we interest the constraint on the 
mass spectrum that is imposed by the experimental results of rare B processes, so we take 
the range of the parameters as: 


-500GeV <B,Bs< 500GeV, 

lO^GeV^ < < 4 x 10®GeV^, 

lOOGeV < mi, m 2 < lOOOGeV, 

-lOOOGeV < /r < lOOOGeV. ( 66 ) 

In the numerical program, the other parameters such as tan/?, tan 6 '„ etc. are given in the 
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figure caption. In Fig. we plot the lightest charged Higgs mass versus 63 (in GeV). Under 
the bound of rare B processes, we hnd the value of Mh+ can vary from lOOGeV to about 
SOOGeV when the parameters vary. Fig. ^ shows the lightest chargino mass varies with the 
parameter 63, with other parameters taken as above. The point we should note is that m-+ 
should be heavy when the tan0„ taken large value (such as tan6*„ ~ 20). When the tan6*„ 
taken small value, the mass of can vary from 30GeV to several hundreds GeV. 

Now, we turn to discuss the branch ratios of rare B processes. From Eq. (B. we hnd 
the parameters B, B 3 enter the mass matrix of charged Higgs (just as in the mass matrices 
of GP-odd Higgs and GP-even Higgs) in forms Bfi and B^e^. Because we are interest in 
relatively light charginos and scalar particles, we take 

-lOOGeV <B < lOOGeV, 

-lOOGeV < H3 < lOOGeV, 

lO^GeV^ < m^i, tuq/, m^/,< 2.5 x lO^GeV^, 

lOOGeV < mi, m 2 < 500GeV, 

-500GeV < p < 500GeV. (67) 

In Fig. 1^ we plot the branching ratio of 6 —>■ s + 7 as a function €3 under some different 
value of tan/3 and tan^^. The dependence on the remnant SUSY parameter such as /i, 
mi, m 2 is represented by the vertical width of the band. We see that positive interference 
with the different sources of SUSY contributions can produce an intensive enhancement 
over the QGD- corrected SM prediction (horizontal solid line) when those parameters are 
assigned suitable values. Sometimes, the supersymmetric contributions dominate over the 
SM contributions resulting in signihcant deviation from the SM prediction. 

Let ns discuss the semileptonic FGNG decay b —>• s + e+e" and 6 —>• numerically. 

In the SM, the QGD-corrected Br{b —>■ s + e’''e“) is about 9 x 10“®. The addition of the 
one-loop contributions where SUSY particles are present modify the prediction up to about 
four times the SM prediction as it can be gathered from Fig. The GLEO collaboration 
has already been searched for inclusive b = e,/i). The results are 
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(68) 


Br{h —> se^e < 5.7 x 10 

Br{h < 5.8 x 10“^ 


By comparing the numerical result with experiment, we hnd that we can not excluded the 
large value of 63. Analogous considerations hold for h (Fig- i|)- The R-breaking 

terms can have an appreciable effect for the b —>■ sJ2 This is the main difference between 
the BRPV model and the usual SUSY model with R-parity. 

In summary, as a simple extension of the MSSM which introduce R-parity violation, the 
R-breaking terms in BRPV model can give an appreciable effect for the rare B-processes. 
From the present experimental bound on those processes, we can get some constraint on the 
mass spectrum in this model under some suitable assumptions. 


Acknowledgment This work was supported in part by the National Natural 
Science Foundation of China and the Grant No. LWLZ-1298 of the Chinese Academy of 
Sciences. 


APPENDIX A: THE DEFINTION OF VARIOUS FUNCTIONS 


We collect in this appendix the various functions that were used in the text. The one- 
variable functions obtained from the penguin diagrams are given as: 


f2\x) = 

ffa{x) = 
= 

fS’ix) = 


X 


1 — X 
1 


1 — X 
1 


Inx, 

{1 + 


X 


X 


{1 + 


1 — X 
1 


Inx}, 


X 


2 In a; 
1 — X 


-|- Inx — 


p(i), . 2a: In a: , a;^lna: 


= 


— Inx -|- 


Ina:}, 

In a: 1 

(1 — xy (1 — x) ’ 

X 


+ 


X 


[l — xY (1 —x)’ 
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c(l) / 
4a I 

p(l) / 


1 .1 


X 


X 


(a; - 1)^2 ^ 

a; — 1 

1 

1 

to 

In 

- ^ 1 - 

X 

X 

In 

(a; - 1)^2 

X — 1 

{x — 1)2 

- 1 rl. 

1 

1 



{x — 1) 2 X — 1 (a: — 1) 


2 Inx], 


ftai.x) = 

fS\x) = 

= 

/m (a^) = 


+ 


X 


+ 


X 


X 


3(a; — 1) 2(a; —1)2 (x — 1)^ (x — 1)^ 


Inx, 


3a; 3a;2 

+ 


3a;2 


2(a;-l) 2(a;-l)2 (a; - 1)3 (a; -1)^ 

1 3a; 3 3a; 

+ -TTTT + 


2(a; — 1) 2(a; — 1)^ (a; — 1)3 [x — 1)' 


In a;, 


In a;, 


1 1 
+ 


1 1 

+ ^- TTT In a;. 


3(a; — 1) 2(a; — 1)^ (a; — 1)3 (a; — 1)^ 

The two- and three-variable functions obtained from penguin and box diagrams are 

1 - * V 

f3l\x,y) = 


(Al) 


fidix,y) = 


x-y 

-Ina 

^a;- 1 

1 

.^2 

X 

2(a;- 

yYx-l 

1 [- 

In 


y -1 


In?/], 


Ina; — 


y 


y-i 
1 


In?/], 


X — y{x — ly a; — 1 


f4e\x,y,z) = 


X 


-{x^ y)], 


3x 




(A2) 


APPENDIX B: THE RELEVANT FEYNMAN RULES IN SUSY MODEL WITH 

BILINEAR R-PARITY VIOLATION 

In this appendix, we give some relevant Feynman rules in the supersymmetric model 
with bilinear R-parity violation that were used in the paper. 

It is hrst convenient to introduce the mixing matrices relative to the scalar-quark sector. 
We denote with the squark current eigenstates (where / = 1, 2, 3 is the generation label 
and Q = U, D), and with ( 5(2 fhe corresponding mass eigenstates of the I-th generation 
(we neglect the generation-mixing of the scalar-quark and scalar-lepton). The 2x2 mixing 
matrices Zqi are dehned by 
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Qi = 

Qr = z^Ql 


(Bl) 


Similar, the 2x2 mixing matrices Z^i are defined as 


= Z^Ei 


l,i Tpl 
El^i 1 

R^* = Zj’iE/ (B2) 

with U, W are the slepton cnrrent eigenstates and El 2 are the corresponding mass eigen- 

, fi, f 2 ) np qnark and down qnark conplings can be 


states. The S'*, (5'j=i^2,3,4 = G , H 
written as (Fig. ^ 


eml 


emt 


-ZHPl + ^- 

V 2Mw sin 6w sin 9^ sin (3 


Z]IPr\Cjj. (B3) 


P~l= -^ 

\/2Mw sin 6 w sin 6 y cos /3 

The conplings of down qnark, np scalar qnark and chargino are given by (Fig. [10|) 

^l( ^ y'li* _|__ i^2i*y+\p I 

sin 6 w \/2Mw sin 9w sin 9^ sin (3 

T 


eml 


fyli* rv— jj I/O* 


\/2Mw sin 9w sin 9^ cos {3 
The conplings of np qnark, down scalar qnark and chargino are given by (Fig. EP 

-Z%,Zy)P,^ + 


(B4) 


au ^ Z^^iZ~^ H _ ^"‘‘d _ j 

sin 9w y/2Mw sin 9w sin 9^ cos (3 

T 


eml 


emt 


yli y+* p -]pi* 


(B5) 


\/2Mw sin 9w sin 9^ sin (3 
The conplings of and Si are (Fig. P^ ) 

~ - ZhZI^*]{p + kY, 

2 sm 9w cos 9w 

where Si = (G“, H~, ti, r 2 ), i = 1, 2, 3, 4. The matrices Z^^ Z^ have been defined as 
above. 


(B6) 


APPENDIX C: THE MIXING BETWEEN r-NEUTRINO AND NEUTRALINOS 


In this appendix, we give the mass matrix of r-nentrino and nentralinos nnder onr 
assnmption. 
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From Eq. (j^) and Eq. (|^, the r neutrino-neutralino mass terms in the Lagrangian are 


C-neutralino = --+ h.C. 


(Cl) 


with = {-iXb, -iX\, - 0 )^ 1 , - 01 ^ 2 , i^tl) and 


2mi 

0 

-\g'vi 

\9'V2 

-\9'v^ 

0 

2m2 

\9Vi 

-\9V2 

\9V3 

-^vi 

\gvi 

0 

-1/i 

0 

yV2 

-\9V2 

-1/i 

0 




0 


0 


The formulae of mixing matrix are: 

- jAb = Zli'x", 

IAa — Xi ) 

/1 0 
i’m — Xi ! 

/ 2 '7^1^ 0 

Xi 1 
= Zn Xi 


and 


= 


Vx?y 


(C2) 


(C3) 


(C4) 


K, = 


X?+l 


(C5) 


y Xi+i j 

with i = 1, 2, 3, 4. Here, we identify the as the lightest mass eigenstate of the mass 
matrix. The matrix Zjsf satishes the following condition: ZJjMnZn = diag{mu^, m^o, m^o, 
m^o, m^o). Similar to the mixing matrices of r-chargino sector, we can assume rriy^, m^o 
{i = 1, 2, 3, 4) positive and m^o > m^o > m^o > m^o > 
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FIGURES 



FIG. 1 . The Feynman diagrams that contribute to 657 and bsZ^ coupling 
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FIG. 4. Considering the experimental results of 6 —> s-l- 7 , 6 ^ se'’"e“, Mfj+ vary with £3 when 

(a)tan0^ = 20, tan/3 = 2; (b)tan0^ = 20, tan/3 = 40; (c)tan0^ = 0.5, tan/3 = 2; (d)tan0„ = 0.5, 

tan/3 = 40. The other parameters are given in Eq. (^). 
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FIG. 5. Considering the experimental results of 6 —> s + 7 , 6 ^ se'^e~, m + vary with £3 when 
(a)tan0^ = 20, tan/3 = 2; (b)tan0.^ = 20, tan/3 = 40; (c)tan0.^ = 0.5, tan/3 = 2; (d)tan0„ = 0.5, 
tan/3 = 40. The other parameters are given in Eq. (^). 
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tan /3 = 40. The solid-lines are the predictions of SM. 
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FIG. 9. The coupling between Si = {G , H , ) and qj 
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FIG. 10. The coupling between fi = {T , K 2 , k^) and qj, Uj 
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FIG. 12. The coupling between and Si = {G 
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